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Abstract. We present a perturbative construction of interacting quantum 
field theories on smooth globally hyperbolic (curved) space-times. We develop 
a purely local version of the Stiickelberg-Bogoliubov-Epstein-Glaser method 
of renormalization by using techniques from microlocal analysis. Relying on 
recent results of Radzikowski, Kohler and the authors about a formulation of 
a local spectrum condition in terms of wave front sets of correlation functions 
of quantum fields on curved space-times, we construct time-ordered operator- 
valued products of Wick polynomials of free fields. They serve as building 
blocks for a local (perturbative) definition of interacting fields. Renormaliza- 
tion in this framework amounts to extensions of expectation values of time- 
ordered products to all points of space-time. The extensions are classified 
according to a microlocal generalization of Steinmann scaling degree corre- 
sponding to the degree of divergence in other renormalization schemes. As a 
result, we prove that the usual perturbative classification of interacting quan- 
tum field theories holds also on curved space-times. Finite renormalizations 
are deferred to a subsequent paper. 

As byproducts, we describe a perturbative construction of local algebras 
of observables, present a new definition of Wick polynomials as operator- va- 
lued distributions on a natural domain, and we find a general method for 
the extension of distributions which were defined on the complement of some 
surface. 
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1. Introduction 

The quest of the existence of a non trivial quantum field theory in four space-time 
dimensions is still without any conclusive result. Nonetheless, physicists are work- 
ing daily, with success, on concrete models which describe very efficiently physics 
at wide energy scales. This description is based on expansions of physical quanti- 
ties like amplitudes of a scattering process in terms of power series of "physical" 
parameters, as coupling constants, masses, charges. The higher order terms of 
these power series are usually ill defined, in a naive approach, but physicists have 
learned soon how to make sense out of them through the procedure now known as 
renormalization |5lj ] . The rigorous extension of this procedure to curved Lorentzian 
space-times will be the main topic of this paper. The question whether the power 
series approximate the corresponding quantities in a full quantum field theory goes 
beyond the scope of this paper and will not be touched. 

The renormalization procedure on Minkowski space-time led to impressive results 
in the case of quantum electrodynamics [ p8L |62| , where observable quantities were 
calculated and agree with high precision with the experimental values . Based 
on this example, a general method of renormalization of interacting fields was found 
and successfully applied to the standard models of elementary particles. 

There is another approach to quantum field theory ("axiomatic quantum field 
theory") which assumes the existence of a class of models satisfying certain first 
principles. Under this assumption several structural properties could be derived 
which therefore hold for every model in this class. To name some, CPT and Spin- 
Statistics Theorems are among the main successes of this line of thought, and nowa- 
days the application of these methods to specific kind of theories, like conformally 
covariant theories on low dimensional spaces, is expected to give new insights. Nev- 
ertheless, all these schemes (see for instance |l^] for a recent survey), either analytic 
[55| or algebraic [^9| , by now seem to have missed the challenge for the concreteness 
needed by, say, particle physicists. 

A notable exception is the rigorous formulation of perturbation theory |34|, 
j53[ |64| ] which may be considered as interpolating between the world of phenomeno- 
logical physics and the mathematical schemes mentioned above. This point of view 
has been pioneered, in particular, by K. Hepp who gave solid foundations to 
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perturbation theory for quantum field theories on Minkowski space-time. This phi- 
losophy proved to be correct for instance in constructive quantum field theory [ p6| 
where rigorous renormalization ideas were used as fundamental inputs. 

One of the aims of the authors is to put forward a formulation of perturbation 
theory which satisfies the needs of axiomatic field theory, much in the sense of 



53 1 , and is at the same time applicable to phenomenology In distinction to earlier 
approaches we give a purely local formulation which is meaningful also on curved 
space-times. 

Our principal result is: 

Main Theorem . All polynomially interacting (scalar) quantum field theories on 
smooth curved globally hyperbolic space-times of dimension four follow the same 
perturbative classification as on Minkowski space-time. 

Before starting the description of our claim we continue the description of the 
interplay between perturbation theory and rigorous methods. 

One of the most puzzling things in physics is that all the attempts to include 
gravity in the renormalization program failed: More recent proposals look for theo- 
ries of a different kind like string theory p^j or its generalizations which are hoped 
to describe all known forces, or Ashtekar program (![]. 

Because of the large difference between the Planck scale (« 10 -33 cm), where 
"quantum gravity" effects are expected to become important, and the scale relevant 
for the Standard Model (3l) (» 10~ 17 cm), a reasonable approximation should be to 
consider gravity as a classical background field and therefore investigate quantum 
field theory in curved space-time. This Ansatz already led to interesting results, 



the most famous being the Hawking radiation of black-holes 32 . A look through 
the literature (see, e.g., |30| ||, H|), however, shows that predominantly free field 
theories were treated on curved physical backgrounds. In fact, most of the papers 
on interacting quantum field theories on curved space-times deal with the (locally) 
Euclidean case and discuss renormalization only for particular Feynman diagrams. 
We are aware of only one attempt to a complete proof of renormalizability, that 
given by Bunch [jl3| for the case of A</> 4 -model. His attempt was however confined 
to the rather special case of real analytic space-times which allow an analytic con- 
tinuation to the (locally) Euclidean situation. It is interesting to note that the main 
technical tool of his paper is a kind of local Fourier transformation and that some 
of his mathematical claims can be justified in the framework of Hadamard "para- 
metrices," both of which belong to the powerful techniques of microlocal analysis 
that we use in this paper. 

The situation is then uncertain for general smooth space-times with the Lore- 
ntzian signature for the metric. Here, to our knowledge, more or less nothing has 
been done. 

Why is the problem of renormalization so difficult on curved space-times? 

Again the precise perspective gained from the rigorous approach is helpful. The 
main problem is the absence of translation invariance which in the rigorous schemes 
plays a decisive role. In general, no global (time-like) Killing vector field exists (no 
energy-momentum operator), so there is no canonical notion of a vacuum state, 
which is a central object in most formulations of quantum field theory; the spectrum 
condition (positivity of the energy-momentum operator) can not be formulated. 
There is no general connection between quantum field theories on Riemannian and 
Lorentzian space-times corresponding to the Osterwalder-Schrader Theorem p7j, 
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and the meaning of path integrals for quantum field theories on curved space-times 
is unclear. As a result, the rigorous frameworks described before cannot simply be 
generalized, and the more formal description based on Euclidean ideas and path 
integrals does not help much. 

On the other hand, physically motivated by Einstein Equivalence Principle p9| . 
a quick look at the possible ultraviolet (short distance) divergences indicates that 
they are of the same nature as on Minkowski space, so no real obstruction for renor- 
malization on curved space-times is visible. Despite of the interest in its own right, 
renormalization on curved space-time might also trigger a conceptual revisitation of 
renormalization theory on Minkowski space in the light of the principle of locality 

m 

To develop perturbation theory in a form which is suitable to extensions to curved 
Lorentzian space-times, we mainly rely on a construction given by Epstein and 
Glaser [ ^p| at the beginning of the seventies and on some improvements suggested 
later by Stora This construction makes precise older ideas of Stiickelberg [[56) 
and Bogoliubov p]. In spite of its elegance it was widely ignored (compare, e.g., 
its neglection in several books on quantum field theory, with the exception of fj8|| ). 
Recently, it was further developed and applied to gauge field theories by Scharf and 
his collaborators |3(j) [] (after earlier work by 0). 

We offer some intuitive explanations of the ideas behind the approach of Stii- 
ckelberg-Bogoliubov-Epstein-Glaser (we refrain from using an acronym for this and 
simply call it the Epstein and Glaser approach). For simplicity we discuss the case 
of flat Minkowskian space-time. 

The basic idea is that, in the asymptotic past and future, the interacting quantum 
fields approach, in a sense to be specified, free fields, i.e. fields satisfying linear 
hyperbolic equations of motion. For free fields there exists a precise construction 
which can be used for a perturbative description of interacting fields. Now, in a 
translationally invariant theory, the interacting fields approach the asymptotic free 
fields only in a rather weak sense (LSZ-asymptotic condition j34|). Moreover, Haag 
Theorem J2i| forbids the construction of interacting fields in the vacuum Hilbcrt 
space of the time-0 free fields. In the Epstein and Glaser scheme, these problems 
are, in a first step, circumvented by choosing interactions which take place only 
in a bounded region of space-time. Then the scattering operator can be defined 
in the interaction picture as the time evolution operator from the past, before the 
interaction was switched on, to the future, after the interaction was switched off. 

A localized interaction here is thought to be a smooth function of time, t, with 
compact support with values in the local operators associated with the free field. 
In the simplest case it is Hi nt (t) = ip(ft) where tp is a free field, i.e. an operator- 
valued distribution on a Hilbcrt space, and where ft{x) — S(x° — t)f(x) for some 
test function /. The S'-matrix is then an operator-valued functional S(f) on the 
test function space. The functional equation for the evolution operator implies 
a factorization property for the S'-matrix if the support of the interaction (as a 
function of time) consists of disjoint intervals. In the case above with the interaction 
being a linear function of the free field we even find the factorization, 

S(h + h) = S{h)S{h) , (i) 



1 Note that the term "finite" in Scharf 's book refers to the fact that in the Epstein and Glaser 
approach (as in the similar BPHZ method) no regularization is necessary. It does not mean that 
the indeterminacy connected with the divergence of naive perturbation theory disappears. 
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whenever there exists some t £ R such that supp(/i) C {x\x° > t} and supp(/2) C 
{x\x° < t}, and where fi and f 2 need not be smooth at the hypersurface x° = 
t. This stronger factorization property is not expected to hold for more singular 
interactions. Instead we require the following consequence of (111), 



to hold for test functions fi, f 2 ,f3, whenever the supports of /i and fa can be sep- 
arated by a Cauchy surface such that supp(/i) lies in the future and supp(/a) in 
the past of this surface. Together with the normalization condition S(0) = 1 (iden- 
tity operator on Hilbcrt space) it implies the first mentioned weaker factorization 
condition in the case f 2 = 0. 

The functional equation (|J) has an interesting property; if S is a functional 
solving it we get other solutions Sf by defining Sf(g) = S(f)~ 1 S(f + g) (the 
relative .S'-matrices) where / is an arbitrary test function. In particular we get 
commutativity in case supp(gi) and supp(<?2) are space-like separated, 



Thus the relative S'-matrices satisfy the locality condition required for local observ- 
ables. They serve as generating functionals for the interacting fields. 

Unfortunately, a construction of S(f) in four dimensions is known only in the 
case of interaction Hamiltonians which are linear or quadratic in the free field, 
but in two dimensions Wrezinski [|33[ proved that, at least in the particular case 
of factorizable /, f{t,x) = g(t)h(x), such a construction is possible for (p 4 . One 
therefore mainly relies on the "infinitesimal" description of the local S-matrix S(g) 
by studying its formal power series S expansion in terms of the "coupling constant" 
g. The connection with the usual formulation may be done via the adiabatic limit, 
i.e. the limit for S(g) when g — > 1 all over space-time is the S'-matrix, or in cases 
where the limit for S(g) does not exist due to infrared divergences, the limit for 
the vacuum expectation values of S g (f), g — ► 1, is the generating functional for the 
time-ordered correlation functions. 

The description given so far emphasizes the fact that the Epstein and Glaser 
method is local in spirit, and it might be a favorite candidate for developing a 
renormalization theory on curved space-times. A closer inspection, however, shows 
that also in this method translation invariance plays an important role, both con- 
ceptually and technically, and it will require a lot of work to replace it by other 
structures. A similar problem was studied by Dosch and Miiller jl6|. These authors 
developed the Epstein and Glaser method on Minkowski space for quantum elec- 
trodynamics with external time independent electromagnetic fields. Their use of 
the Hadamard parametrices for the Dirac operator is already much in the spirit of a 
local formulation of perturbation theory; by the assumption of time independence 
of the external fields, however, time translation invariance still plays a crucial role 
in their approach. As a matter of fact, it will turn out that microlocal analysis [ |35| 
is ideally suited to carry through the program where in particular the concept of 
the wave front set proves to be extremely useful. We note, en passant, that other 
reseachers |l5|, |36|, |37| had previously used these tools in quantum field theory and 
that more recently Verch j58| has developped a generalization of the concept of 
wave front sets which can be applied in algebraic quantum field theory. 

This paper is an extended version of a previous one Q where we sketched the 
main ideas. Here on we give all the necessary mathematical details. 




(2) 



Sf(gi+92) = S f (g 2 )Sf(g 1 ) = S f (g 1 )S f (g 2 ) ■ 



(3) 
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The paper is organized as follows: After this introduction, Section 2 provides 
some useful grounding concepts and fixes the notations. Moreover, we present a new 
construction of Wick polynomials which may be of independent interest. In Section 
3 we state the first principles by which we build up the perturbative method on 
smooth curved globally hyperbolic space-times. The most important change w.r.t. 
the Epstein and Glaser method is a characterization of the singularity structure of 
the time-ordered numerical distributions replacing translation covariance. In the 
course of this part we show a local version of the so called "Theorem 0" of Epstein 
and Glaser which provides the necessary mathematical properties of the building 
blocks of the construction. In Section 4 we start the inductive procedure which aims 
at constructing the time-ordered functions up to the small diagonal of the product 
manifold M n , where all "dangerous" singularities are located. Sections 5 and 6, 
have a more mathematical flavour; we introduce the concept of scaling degree at 
a point, following essentially Steinmann |^3| , and its generalizations in terms of 
microlocal analysis. The main aim of this section is the description of the extension 
to all space of distributions defined on the complement of a submanifold. These 
tools are needed for the classification and implementation of renormalizability. The 
next, Section 7, contains the end of the inductive procedure by which we prove the 
theories with polynomial interactions to follow the same perturbative classification 
as on Minkowski space-time. 

We emphasize that the method of defining the local S'-matrix joins perturbation 
theory with the more abstract algebraic formulation of quantum field theory. In 
fact, we are able to define a unique family (net, precosheaf etc.) of ^-algebras of 
observables on globally hyperbolic space-times via the idea of the local relative S- 
matrices. Section 8 describes this construction which seems to be widely unknown, 
in spite of the fact that it may already be found, in a preliminary form, in j52|. 
This Section partly justify the rather abstract starting point of Section 2. 

An outlook, Section 9, concludes the paper. Finally, we stress that the proce- 
dure works for general field theories but for simplicity we stick to the notationally 
easiest case of a single scalar (massive) field theory with self interactions without 
derivatives. 

2. General Theory of Quantized Fields and Microlocal Analysis 

In order to fix our notations we recall some basic geometrical concepts. Further 
details may be found in some books on general relativity and Lorentzian geometry 
(see, for instance, |5^] and Q). We shall work on a space-time (M, g), where by 
this we mean that M is a connected, Hausdorff, boundary less topological space of 
pure dimension d > 2 which (i) is paracompact, (ii) is equipped with a smooth 
structure, (Hi) is endowed with a Lorentzian metric g a b, i-e. a smooth 2-cotensor 
of signature (1, d — 1), i.e. (+, — , • • • , — ) and (iv) is oriented and time oriented. 
Given the metric we have a canonically associated derivative, namely the Levi- 
Civita derivative denoted by V and an associated curvature tensor R^ cd , with R 
the scalar curvature. The notion of totally geodesic submanifold, i.e., that one for 
which all tangential geodesies stay on the submanifold, is used in Section 6. 

Some words on notations. Sometimes we write a zero section of a vector bundle 
® as {0} some other times to precise that it belongs to that bundle we shall write 
Z(H). However, in order to avoid any abuse, we use the notation 23 to denote the 
bundle deprived from its zero section, i.e., 23\Z(23). We shall also use the notation 
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M n , whenever we treat the n-th order cartesian product of a manifold M, and by 
Ak, where K C {1, . . . ,n}, we mean the (smooth, closed) submanifold of M n for 
which any of its points (pi, . . . ,p n ) are such that pk ± — Pk 2 f° r an y P a i r &l 7^ ^2 in 
if. 

The causality principle plays a crucial role in our construction. Therefore we 
restrict our space-times to be globally hyperbolic. This means that M is homcomor- 
phic toKxS where £ is a (d — l)-dimensional topological submanifold of M and 
for each t € M, {i} x S is a (spacelike) Cauchy surface. A Cauchy surface is a subset 
of M which every inextendible non spacelike curve intersects exactly once. Given 
a subset S of M we define the causal future/past sets J ± (S') as the subsets of M 
which consist of all points p 6 M for which there exists some point s € S connected 
to p by a non space-like future/past directed curve. If M is globally hyperbolic, 
the set J + (p) n J~(q) is compact for any pair p,q & M. Finally, if p e M then the 
induced metric on tangent space T p M and cotangent space T*M are Minkowskian, 
and we define the future/past light-cones V± over these spaces (based on p) in the 
usual way. 

Quantum field theories on more general spaces pose consistency problems (see, 
e.g., Hawking's "Chronology Protection Conjecture" Q or the divergence of the 
energy momentum tensor at the Cauchy horizon observed in |40|). We remark, 
however, that since our constructions will be purely local, one can as well consider 
a globally hyperbolic submanifold of any Lorentzian space-time. 

In many concrete cases, exact solutions of the Einstein equation, like Minkowski, 
de Sitter, Schwartzschild are real analytic. In these cases some of our results might 
be sharpened by working with the analytic version of microlocal analysis H] . In 
this respect, we should mention some recent results of Bros, Epstein and Moschclla 
for a Carding- Wightman-like description of quantum field theories on de Sitter 
space-time where analytic function techniques play a major role. 

2.1. Wave front sets and Hadamard states for free fields. For the (massive) 
free field ip satisfying the (generalized) Klein-Gordon equation of motion 

(D g + to 2 + nR)<p = , (4) 

where D 9 is the d'Alembertian (or Laplace-Beltrami) operator w.r.t. the Lorentzian 
metric g, m > and k £ K, one may associate an algebra of observables defined 
in the following way: Let E ret resp. E a d v be the retarded resp. advanced Green 
functions of the Klein-Gordon operator which are uniquely defined on globally hy- 
perbolic space-times, and let E = E ret —E ac iy Then we consider the unital ^-algebra 
21 which is generated by the symbols <p(f), f G D(M) (space of complex-valued 
smooth and compactly supported functions), with the following relations: 

1. The map / i— » <p(f) is linear, 

2. ¥>(/)* = <p(J), 

3. Mf), ip(g)} = iE(f 5 )1, V/, g e D(M), 

4. <p((D g + m 2 + K R)f) = 0, V/ S D(M), 

where the symbol [</?(/), f(g)] stands for <p(f)<p(g)—<p(g)<p{f) and / means complex 
conjugation. A state is, by definition, a linear functional to on 21 (the expectation 
value) which is positive (i.e. oj(a*a) > 0) and normalized (u)(l) = 1). It is uniquely 
determined by a sequence of multilinear functionals uo ni n — 0, 1, . . . (the n-point 
functions) on the test function space D(M), 

U) n (f 1} ... ,f n ) = u((p(fi) • • • (p(f n )) ■ (5) 
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We only consider states whose n-point functions are distributions and restrict fur- 
thermore our attention to the states called quasi-free, namely, those states whose 
only non trivial n-point functions have n even and are generated in terms of the 
2-point functions (see, e.g. p9j). Among them a distinguished class is formed by 



the so called Hadamard states (see, e.g. |14| |4l)). They are thought to be the 
appropriate analogue of the concept of the vacuum which has no direct counterpart 
on generic space times. In fact, they are quasifree states whose 2-point functions 
have a prescribed short-distance behaviour which is partially motivated by the fact 
that it allows the definition of the expectation value of the energy-momentum ten- 
sor (see, e.g. |6Cj). As first observed by Radzikowski |l9j the 2-point functions of 
Hadamard states can be characterized in terms of their wave front set. 

To discuss this characterization we need to enter into the realm of microlocal 
analysis. We give some motivations to the basic notions of wave front sets and 
present those basic results which are used throughout the paper. We leave the 
reader the task to look further into the large existing literature |3jJ. Physicists 
might, for concreteness, start from the well- written short exposition of Junker in 
|39| , where they can find definitions and results about pseudodifferential operators, 
which wc hold as known. 

We shall denote by £(R n ) the space of complex- valued smooth functions and by 
£'(R n ) its dual space, i.e., the space of compactly supported distributions. 

It is a standard result in distribution theory that u € £'(R n ) is a smooth function 
iff its Fourier transform u decays rapidly in Fourier dual space R„, i.e. for any 
integer N there exists a constant Cn such that \u(k)\ < Cjv(1 + \k\)~ N for all 
k G R„, where R„\{0} = R„. In case u is not smooth the Fourier transform may 
still rapidly decay in certain directions. We may describe this set of directions by 
an open cone in R„ and define S(u) as its complement in R„. It is easy to see 
that E(<fm) C S(u) when u G £'(R") and (j> S D(W l ). This property suggests a 
strategy for the general case in which u is not of compact support. So, considering 
u G D'(R") and a point x G R n in the support of u, supp(u) C O, O open 
subset of R™, we first localize u via multiplication with some </> G D(R n ) such that 
4>{x) ^ and then consider the Fourier transform of 0u, now a distribution of 
compact support. We then define the set T, x (u) = n^E^u) where the intersection 
is taken w.r.t. all smooth functions of compact support (j) such that (j)(x) ^ 0. 
This may be called the set of singular directions of u over x. It is empty whenever 
x singsupp(u). 

Hence, finally, we define the {smooth) wave front set for u G D'(W ) as WF(u) = 
{(x, k) G R" x R„ | k G S x (w)}. This set is readily seen to be closed and conic, 
where the last means that if k G ^ x (u) so do any Xk for all A > 0. 

It is now crucial that the notion of the wave front set can be lifted to any smooth 
manifold M where it is invariantly defined as a subset of the cotangent bundle T*M. 
This covariance under coordinate transformations is what gives to the definition its 
real technical power. 

Among the results which will be important for us we mention that derivatives 
do not enlarge the wave front set of a distribution, i.e. WF(du) G WF(u), and the 
following criterion called Hormander criterion for multiplication of distributions: 
Ml. Product. Picking two distributions ui,U2 G D'(M), the pointwise product 
u\Ui exists as a bona-fide distribution whenever WF(ui) + WF(tt 2 ) does not 
intersect the zero section Z(T*M), i.e., if for all covectors ki G WF(wj), 
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i = 1,2, based over the same point one finds that k\ + k 2 7^ 0. Moreover, if 
WF(u,) cr„i = 1, 2, then WF( Ul u 2 ) C T 1 U L 2 U (T 1 + T 2 ). 
We shall also refer frequently to a certain continuity property in microlocal anal- 
ysis which in the body of the paper is sometimes called "Hormander (pseudo) 
topology." It has to do with the notion of convergent sequences which respect also 
wave front set properties: 

M2. Continuity. Let T>' r (M) = {v £ D'(M) | WF(w) C T}, where T is a closed 
conic set in T*M. A sequence {iti}i G N C T>' r (M) converges to u £ D r (M) 
in the sense of the Hormander (pseudo) topology whenever the following two 
properties hold true: 

(a) m — > u weakly* (i.e. in D'(M)), 

(b) for any properly supported pseudodifferential operator A such that /J,supp(^4)n 
T = 0, we have that Aui — > Au in the sense of £(M). (/isupp(A) is the projec- 
tion onto the second component of the wave front set of the Schwartz kernel 
of A.) 

A last property is connected with the sequential continuity, in the sense of M2, 
of the operation of restriction of a distribution to a submanifold: 
M3. Trace. Let N C M denote a submanifold, and let u £ D'(M). Then u can 
be restricted to the submanifold 3\f whenever WF (it) does not intersect the 
conormal bundle N*N of N. Moreover, if WF(w) C L, with L a closed conic 
set such that T D N*N = 0, then the operator of restriction (trace) 7 can be 
lifted as a sequentially continuous operator, in the sense of M2, from T>' r (M) 
to D'(N). 

For later purpose, it is convenient to have a coordinate dependent formulation 
of M2(b) by using Fourier transforms. Namely, let xq £ M and let V be an open 
conical neighbourhood of T Xo , where the last denotes a set of covectors associated 
to the point xq. Choose a chart (<p,U) at xq such that T x C V for all x S U. 
Let x S D(U) with x( x o) 0- Then the Fourier transform of \ u , u S T>' r (M), is 
strongly decreasing in the complement of V, and 

sap\0&-xu)(k)\(l + \k\) N ^0, (6) 

k^V 

for all N £ N if Ui — > m in Dp(M). If, on the contrary, the above convergence holds 
true for all choices of xo, V, (ip, U) and x> we obtain M2(b). 

After this disgression into microlocal analysis we briefly describe Radzikowski's 
characterization of Hadamard states . The idea is to use wave front sets for a 
formulation of a spectral condition. The antisymmetric part of the 2-point function 
is the commutator function E. Its wave front set is 

WF(£) = {(a?, k; x', —k') £ f*M 2 \ (x, k) - (x', k')} . (7) 

Here the equivalence relation ~ means that there exists a null geodesic from x to x' 
such that k is coparallel to the tangent vector of the geodesic and kl is its parallel 
transport from x\ to x 2 - For coinciding points, the relation is defined as consisting 
of the degenerate (i.e., only one point) geodesic at x = x' which has covector k still 
along the boundary of the light-cone and k' = k. We remark, for later purpose, 
that since only light-like covectors are present, one can restrict E, and, whenever 
local coordinates are chosen, its derivative w.r.t. time E, to any spacelike Cauchy 
hypersurface. 
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As a result of |49| [43| , the 2-point function of a Hadamard state has a wave front 
set which is just the positive frequency part of WF(E), 

WF(w 2 ) = {{x, k; x', -k') G WF(£) | k G V+} . (8) 

Since (^|) restricts the singular support of ^(^l, ^2) to points x\ and X2 which 
are null related, 0J2 is smooth for all other points. The smoothness for space-like 
related points is known to be true for quantum field theories on Minkowski space 
satisfying the spectrum condition by the Bargmann-Hall-Wightman Theorem |]55| . 
For time-like related points, however, a similar general prediction on the smoothness 
does not exist. 

Another deep result from Radzikowski jl9| shows that the Duistermaat-Horman- 
der [p~7| distinguished parametrices for the Klein-Gordon equation are nothing else 
than the (Stuckelberg-)Feynman-anti- Feynman propagators (up to C°°) for quasi- 
free Hadamard states. We recall that the time-ordered 2-point function Ep arising 
from u>2 is given by 

iE F (xi,x 2 ) = uj 2 (xi,x 2 ) + E ret (x 1 ,x 2 ) ■ 



Its wave front set [49 



is 

WF(E F ) = OUD , 

where the off-diagonal piece is given by, 

O = {{x,k;x',~k') G T*M 2 \ (x,k) ~ (x',k'),x^ x',keV± if x G J^z')} , 

and the diagonal one by, 

D = {{x, fc; x, -k) G T*M 2 \x&M,ke f*M} . 

Now, one can see why in naive perturbation theory we may find divergences. 
Indeed, the perturbative expansion in terms of Feynman graphs in position space 
leads to pointwise products of Feynman propagators. But these products do not 
satisfy Hormander criterion for multiplication of distributions since covectors based 
on the diagonal piece D can add up to zero. 

2.2. A new construction of Wick polynomials. In a previous paper || we 
constructed Wick polynomials as operator- valued distributions. We considered a 
fixed Hadamard state to and the induced GNS representation (Ti.^ , ir u , il^ ) for the 
*-algebra 21 and found the Wick polynomials as operator- valued distributions on the 
dense cyclic domain generated by O w . We recall that by a GNS triple {Ti. w , n u , ) 
we mean a complex Hilbert space Ti.^,, a representation 7r w of 21 by unbounded 
operators on 7i w , and finally by fl^ £ Hu> the cyclic vector representing the state 
uj for which one has the connection equation 

u(A) = (n^TT^A)^) , VAG21. 

The dependence of the construction of Wick polynomials on the choice of the 
Hadamard state led to two problems: The first one is due to the convention that 
the expectation value of a Wick polynomial vanishes in the chosen Hadamard state. 
Other choices lead to a finite redefinition, a problem well known from the definition 
of the expectation value of the energy momentum tensor ]6l[ | . Since we shall not 
discuss finite renormalizations in this paper, we do not treat this problem at the 
moment. The other problem is of a more technical nature: The smeared Wick 
polynomials are unbounded operators. We know from the work of Verch |57]] that, 
locally, i.e. in bounded regions of space time, different Hadamard states lead to 
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equivalent representations. But this theorem does not guarantee that the domains 
of definition for different choices of the cyclic vector coincide. We therefore give 
here a new definition which depends only on the representation but not on the cyclic 
vector. Its restriction to the cyclic subspaces coincides with the previous definition. 

It is well known that the operators <p(/) (now, representatives under tt u of the 
abstract elements of the *-algebra 21 in subSection 2.1) for a real valued test function 
are essentially self adjoint on the cyclic domain generated by n w , and that the Weyl 
operators W(f) = exp(iip(f)**) (where now the *-operation denotes the Hilbcrt 
space adjoint) satisfy the Weyl relation W(f)W(g) = exp(—%E(f ig) g))W{f + g). 
The expectation value in the given Hadamard state is 

w(W(/))=exp(-iwa(/,/)). (9) 

Let :W(f): = exp(jW 2 (/, f))W(f), and define for * 6 H u the vector-valued func- 
tion *(/) = :W(f): * . 

Definition 2.1. We say that ^f(f) is infinitely often differentiable at f = if there 
exists for every integer n > a symmetrical vector-valued distribution S 11 ^ /5f n 
on D(M n ), and continuous seminomas p n on the test function space T>(M) with 
Pn+i > Pn, such that 

(a) ifp (h) = then V(h) = *(0) , 

(b) if h — > with p n {h) ^ then 



1=0 ' J 



PnQl)- n — » 



where \\ ■ \\ stands for the Hilbert space norm in T~L U . 
The kernel of the functional derivative can be written, 

— j— — -^=i n :ip( Xl )---ip(x n ):* . (10) 

ojyxi) ■ ■■of{x n ) 

The right hand side of Eqn. ( p^| ) defines what is called a Wick monomial. 

We want to find those vectors on which the Wick monomials can be restricted 
to partial diagonals. In view of the criterion M3 for the restriction of distributions, 
we define as the microlocal domain of smoothness the following set: 

V = |\I> e Huj | *(/) is infinitely often differentiable at / = 0, and for every 

ri 6 N the wavcfront set of — — is contained in the set 

Sfn _ , (ID 

{(x 1 ,k 1 ;...;x n ,k n ) | h G V-, i = l,...,n} 

The vector-valued distributions (|l0|) with $ e P can be restricted to all partial 
diagonals, and give all possible Wick polynomials. Moreover, according to Ml, 
they may also be multiplied by distributions whose wavefront sets do not contain 
elements (x±, fci; . . . ; x n , k n ) € T*M n , ki £ V+, i = l,...,n. The domain T> is 
invariant under application of Weyl operators and smeared Wick polynomials. 

We state a crucial property: 

Lemma 2.2. Let $ € Ti.^ induce some quasi- free Hadamard state u>'. Then $eP. 
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Proof. The main point rests on the validity of Leibniz rule. Indeed, we can write, 

$(/) = exp(±(a, 2 (/, /) - Lj' 2 (f, /)) exp(i^(/, /) + , 
and differentiate w.r.t. /. The general n-th order derivative gives 

g>*)= E x(h,h)^g^(h^n h (12) 

^ jc{i,...,n},|i'| even J 

where, <&'(/) = exp(iw2(/, /) + i(p(f)**)<&, x — ^2 — ^2 1S a smooth function on 
M 2 and a solution of the Klein-Gordon equation in b oth entries. 

Now, (f>'(/i) satisfies the estimate in Definition 2A with the seminorms p' n {h) = 
y^uj' 2 (h, h) for all n, and the numerical prefactor with the seminorms p n {h) — 
y/ (u>2 + w' 2 ) (h, h), hence for the whole expression we may also use the seminorms 
p n . Actually, as was shown by Verch in f57| , there exist two positive constants A 
and B, such that 

Aua(f,g) < <4(f,9) < Bu 2 (f,9) , T)(M) , 

hence all these seminorms are equivalent. We conclude that $(/) is infinitely dif- 
ferentiate at / = 0. 

The wave front sets of the functional derivatives of <&(/) and <&'(/) coincide, 
since \ is smooth. Using the formula 

2 

= n!w£(M) n , (13) 



■n v 



Sf 

and the information on the wave front set of Hadamard states, we find 
WF T7T ={(x 1 ,k 1 ;...;x n ,k n )eT*M n \ h e dV_, i = 1, . . . 



Sf- 



(14) 



so $ G V. 



We use the so called (generalized) Wick expansion formula which is the basic 
combinatorial formula for perturbation theory. Let us denote by £ any Wick poly- 
nomial on D C Ticj, and we define the "derivative" of a Wick polynomial with 
respect to cp to be dH/dtp. We can characterize it by the following result: 

Lemma 2.3. There is a unique Wick polynomial dL/dip which satisfies the equa- 
tion, 

[Jl(x),v(y)} = —(x)E(x,y), (15) 
in the sense of operator-valued distributions on T>. 

Proof. By linearity it is sufficient to prove it for any Wick monomial e.g. :ip n (x) :. 
It is obvious, that n:ip n ~ 1 (x) : satisfies (filf), hence we only need to proof that if A 
is any Wick polynomial for which A(x)E(x, y) — this means that A = 0. But this 
follows from the fact that for every x 6 M we can find some test function / such 
that the (smooth) solution E(x, f ) of the Klein-Gordon equation does not vanish 
at x. ■ 

Now, let us consider, for any Wick polynomial £, the fields £W = d^H/difP, 
j S N, which, by induction, are uniquely defined according to the previous Lemma. 
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Theorem 2.4 (Generalized Wick expansion Theorem). Let k = 1, ...,n, be 
Wick polynomials. The following relation holds: 



x 




where the summations over the jk 's go from to the order of the corresponding L)~- 

For a proof see p3| or just use the previous notion of differentiability and apply 
induction. Note that the products in the Theorem above exist because the wave 
front sets of their expectation values satisfy Hormander criterion Ml due to the 
convexity of the forward light cone. 

The wave front sets for the Wightman distributions of Wick polynomials may be 
larger than those of the corresponding distributions for the free field ip. Consider 
as an example the 2-point Wightman function for the Wick monomial :ip 2 (x) :, i.e. 

:ip 2 (xi) ::(y9 2 (x 2 ) :f2 w ). According to the Theorem, this is equal to 2w 2 (^i, X2) 2 ■ 
This product exists according to the Hormander criterion, and its wave front set is 
contained in (WF(w 2 ) + WF(w 2 )) U WF(w 2 ) = Fi- The set F 2 will be instrumental 
for some results below. Now, WF(w 2 ) + WF(w2) contains directions which lie inside 
the light-cone as it is clear by adding up two covectors k\ + &2 for points on the 
diagonal. One thus sees how already the smallest possible non-linearity may give 
rise to additional singular directions w.r.t. those already present in the wave front 
sets of the Wightman functions for the original field p. Another important remark is 
that F 2 is an involutive closed cone, i.e. is a closed cone which is stable under sums, 
and, as a straightforward result, it gives that (ui2(x, y)) n still has WFf^J) CF 2 - 

The general structure for multi-point expectation values of Wick polynomials 
can be found as follows. For a more compact notation some definitions from graph 
theory are used: Let Q n denote the set of all finite nonoriented graphs with vertices 
V = {1, . . . , n} and let E G denote the set of edges of a given graph G. Moreover, 
for any vertex i £ V we denote by Ef the subset of edges which belong to the 
vertex i, possibly empty, and by \E G \ their number and similarly by E G the subset 
of edges connecting points i and j, with the obvious relation \Ef\ = J^j l-Ey |- For 
any edge e£E connecting points i and j we use the "source and range" notation, 
i.e. i = s(e) and j = r(e), whenever i < j. 

It is sufficient, by linearity, to restrict ourselves to the treatment of products 
of Wick monomials. Indeed, let us denote by a;™ 1 '"''" 1 ™ the expectational value, 
w.r.t. the GNS-vector Q u for a quasi-free Hadamard state u>, of the product of 
Wick monomials :p mi (xi) : ■ ■ ■ :p mn (x n ) : , and define as G n (m\, . . . ,m n ) the set of 
all graphs G for which all vertices j with irij edges are saturated, i.e. \E? \ = nij. 
Moreover, following we call a triple (£,7, k) an immersion of any graph G C Q n 
into the manifold M whenever, (a) x : V — * M is a map from all vertices i of 
G to points Xi of M; (b) 7 maps edges e e E G to null geodesies 7 e connecting 
points £ s ( e ) and x r ( e y, (c) k maps edges e S E G to future directed covector fields 
kj = k e which are coparallel to the tangent vector 7 e of the null geodesic. Hence, 
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generalizing the set Fi above by, 

F n = l{xi, fci; . .. ; x n , k n ) g t*M n \ 3 G g Q n (in±, . . . , m„) and an 
immersion (a;, 7, A:) of G such that 

kl = X! M^i) - ^ fe e (a;i)| , (17) 

e£E G e£E G 
s(e) — i r(e)—i 

we have: 

Proposition 2.5. TTie wave front set of a;™ 1 '— > m ™ is geometrically bounded as 

WFK" 11 m ») C F„ . (18) 

Proof. Let us define for notational purpose C7„ = Q n (mi, . . . , m n ). As follows from 



Theorem 2.4 the expectation value is, 

w™ 1 '---' m "(xi,...,x„) = ^ w^(xi, . . . ,x n ) = ^2 II w 2(a; s ( e ),a; r(e) ) . 

cee„ Gee„ e 6G (19) 

Considering one graph G in this sum we see that, in explicit form, 

e£G 

= W 2 {xi,X2)\ E °*\w2{xi,Xs)\ E °*\ ■■■W 2 {Xn-l,X n )\ E °-^. 

In the last equality the 2-point distributions should be understood as distributions 
on the product manifold M n . Hence their wave front set is given, when i < j, by, 

\E G -\ 

WF(w 2 * ,j ) c {(2:1,0; ...;Xi, kij; . . . \Xj, -k itj ; . . . ;x n ,0)\(x 2 , fa,j-,Xj, -k it j) g F 2}- 

It is straightforward to see from the last expression that the claim of the Propo- 
sition is correct. Indeed a general covector ki will have the following expression, 

h = -ki,i h-i,i + ki.i+i H 1- ki,n > (20) 

where some of the ki tTa may be zero. Now, as follows from Eqn. (^), to any edge 
e which consists of pair of joined vertices i,j in the graph G there exist points on 
the manifold Xj , Xj and a null geodesic j e connecting them together with a future 
directed covector field k e which is coparallel to the tangent vector of the geodesic 
and is such that, in agreement with Eqn. (p0[), 



fcj — S k e ^x{) S k e (xi) . 

e£E G e£E G 
s{e)—i r(e)— i 

Since this applies equally well to all graphs, and since the wave front set of sums of 
distributions is bounded by the union of the wave front sets, we get the thesis. I 

3. On a Local Formulation of Perturbation Theory 



We have recalled in the introduction the main ideas of the Epstein and Glaser 
formulation of perturbation theory. Here we give the details of our generalization. 
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We start from the Gell-Mann and Low formula for the 5-matrix for quantum 
theories on four dimensional Minkowski space-time M; this means adopting the 
following formal expression, 

5 A =T(e iX ^ liint( - x ' >d4x ) , 

where T denotes the notion of time ordering, Hint, the interaction Lagrangian, is 
some local field and A is the strength of the interaction. Developing in Taylor series 
w.r.t. A gives 

S\ = V" — n— / • •• / T(£ Jint (x 1 )L int (x 2 ) ■ ■ ■ &int(xk))T\d' i x 



% • 



Hence, the perturbative solution to scattering theory is reduced to quadratures once 
one finds the general solution of the time ordering operation inside the integral. For 
noncoinciding points the solution is given by the following expression, 

T(&int(xi) ■ ■ ■ £ji n t\ x n)) = 

6( x n(l) ~~ x ir(2)) ' ' ' @( x n(n-l) — x tt (n) )^int (^(1) ) ' ' ' ^int{ x ir{n)) > 

7re3>i,...,„ ( 21 ) 

where Ti,...^ is the set of all permutations of the index set {1, . . . , n} and 6 is the 
Heaviside step function, 



6(x) 



1, ifx°>0, 
0. otherwise, 



where a; denotes the time component of the points in M. As well known, this 
expression leads to the description of scattering processes by Feynman graphs p9[ . 
Due to local commutativity of the Lagrangian the singularities of the Heaviside step 
function at coinciding times are harmless, as long as all points are different. Unfor- 
tunately, this is no longer true for coincident points, since Hint is an operator-valued 
distribution which cannot, in general, be multiplied by a discontinuous function; if 
one tries to define the products by convolutions in momentum space, this leads, in 
a naive approach, to the occurrence of ultraviolet divergences. 

Several procedures have been found to cope with these singularities. But typi- 
cally they are nonlocal and are therefore not immediately generalizable to the case 
of Lorcntzian curved backgrounds. Better is the situation in Euclidean field the- 
ory (see e.g. p5[ for a generalization of dimensional renormalization to the curved 
case) . 

Let us consider, as possible interaction Lagrangians L, Wick polynomials of the 
free field <p. From Eqns. (^lj) and (|l^) we find 



:ip 3l (xi) ■ ■ -^"(xn): 



jl! ■ ■ -jn ] - 



(22) 



where, however, the expectation value in the right-hand side is a priori not defined 
all over M™, but only over M n \E where S is the union of all diagonals in M". So, 
the main problem is to give a mathematical meaning to this formula on all points. 
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3.1. Formulation of the local S- matrix. The general starting idea, due to Bo- 
goliubov, is to consider the usual Gell-Mann and Low formulation of the ^-matrix 
but supplemented by the hypothesis necessary to the implementation of the causal- 
ity principle. In one stroke one finds also the solution to the problem of the correct 
treatment of the operator-valued distributions. Choosing as the interaction La- 
grangian L;„i(i) = L{x)r}{x), a Wick polynomial L multiplied by a space-time 
function of compact support r\ G D(M) (considered as a generalized "coupling 
constant" ) , we define the local S- matrix S\ (rj) as a formal power series (see, for 
instance, ||) in the coupling strength A, 

00 (iX) n f 

where dfi is the natural invariant volume measure on the globally hyperbolic space- 
time (M,g), and 1 is the Hilbert space identity operator. One can enlarge the 
definition e.g. by using a more general Lagrangian, 

l 

A£;„ t = Vk&k , (24) 
fe=l 

with Wick polynomials Lk and associating to each a different "coupling constant" 
rjk G D(M), where the additional "Lagrangians" are defined as terms like cur- 
rents, external fields etc., including in particular all derivatives of the basic interac- 
tion Lagrangian. Using this extended Lagrangian we may replace the time-ordered 
operator in Eqn. ( p3| ) by, 

where the summation over the fe's go from 1 to I, the number of the terms in the 
extended Lagrangian. 

We remark that, eventually, the test function(s) 77 should be sent to a fixed value 
over all space-time. This procedure, known as adiabatic limit, amounts to treat the 
infrared nature of the theory. Some studies of this limit in the case of Minkowski 
space-time have been performed by Epstein and Glaser themselves [^lj. It is not 
clear how to generalize their study to curved spaces. It is therefore gratifying 
that all local properties of the theory are already obtained via the construction 
of the local S-matrices, and this point of view might also be useful in cases (like 
non abelian gauge theories) where due to infrared problems the S'-matrix in the 
adiabatic limit does not properly exist. 

3.2. Denning properties. Our main goal is the inductive construction of the 
time-ordered products of Wick polynomials, 

T(Li(xi) ■ ■ • . 

Following Epstein and Glaser we require the following properties: 

PI. Well-posedness. The symbols T(Li(xi) ■■■ L n (x n )) are well defined ope- 
rator-valued distributions on the GNS-Hilbert space H u i.e. (multilinear, 
strongly continuous) maps D(M n ) — * End(2?) where T> C 7i u is the dense 
subspace (microlocal domain of smoothness) defined in (|Tl|). 
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P2. Symmetry. Any time-ordered product T(Li(xi) ■ ■ ■ L n (x n )) is symmet- 
ric under permutations of indices, i.e. the action of the permutation group 
< ^{x,...,n} 01 the index set {1, . . . ,n} gives, 

for any n £ IP{i ,„}, in the sense of distributions. 

This symmetry property corresponds to the fact that the time-ordered prod- 
ucts are functional derivatives of the local iS-matrix. More crucial is the following 
causality property, which follows from the Eqn. (^); 

P3. Causality. Consider any set of points (xx, . . . , x n ) G M n and any full parti- 
tion of the set {1, . . . , n} into two non empty subsets I and I c such that no 
point Xi with i E I is in the past of the points Xj with j 6 I c , i.e. Xi (£ J~(xj) 
for any i G I and j € I c ■ Then the time-ordered distributions are required to 
satisfy the following factorization property 

T(H 1 (x 1 )---H n {x n ))=T(l[£ Ji (x i )) T(J] Z,(x 3 )) . 

iel j£l c 

In the Epstein and Glaser scheme on Minkowski space-time one requires, in 
addition, translation covariance of the time-ordered products. If the free field is 
among the possible terms in the Lagrangian one can show that the time-ordered 
products are sums of pointwise products of Wick polynomials with translation 
invariant numerical distributions. (Such products exist due to Theorem of Epstein 
and Glaser, an easy proof of which follows from our microlocal characterization of 
the domain of definition of Wick products.) Moreover, the condition, 

[T(£i(xi) • • • £„(£„)), ip{y)\ = y2 T(£i(xi) ^i x i) ' ' ' &n(x n ))E(xi,y) , 

i=i d(p (25) 

fixes the coefficients to be vacuum expectation values of time-ordered products 
of those Wick polynomials which are of lower order w.r.t. the chosen interacting 
Lagrangian (from now on, we shall call them su&-Wick polynomials), hence the 
problem is reduced to a problem for numerical distributions. Unfortunately, in the 
case of a curved space-time, we have not yet determined the class of fields which are 
relatively local to the scalar free field, i.e., what is known in the literature as the 
Borchers' class [^2|. One also needs a replacement of the condition of translation 
covariance. Our idea is to impose a condition on the time-ordered distributions 
which in a sense employs both ideas of invariance and spectrality crucial in the 
Minkowskian case. Since, as emphasized in the previous section, spectrality for 
us means wave front sets properties we now look for a condition which fixes the 
properties of the singularities of the time-ordered distributions. We use the graph 
theoretic definitions of Section 2. 

P4. Spectrality. For the expectation value t n £ D'(M"), n > 2, of any time- 
ordered product it holds, 



wf(*„) c r 
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where 

C = |(a;i,fci;...;a; n ,A: n ) G f*M n \ 3 a graph G e G n and 

an immersion (x, 7, fc) of G in which k e is future directed 
whenever x a r e ) ^ J — (a; r ( e )) an d such that, 

m:s(m)—i n:r{n)—i 

This may be motivated by the fact that, for non coinciding points, £ n can be 
expressed in terms of the usual Feynman graphs and for the set of coinciding points 
we have an infinitesimal remnant of translation invariance, since all covectors at 
coinciding points sum up to zero. 

We can now formulate a microlocal version of Theorem of Epstein and Glaser. 

Theorem 3.1 (Microlocal Theorem 0). -Z/P4 holds for t n then 

t n (x\, . . .,x n ):(p ll (xi) ■ ■ ■cp ln (x n ): , 

is a well defined operator-valued distribution, for any n and any choice of indices 
l±, . . . , l n , on the dense invariant domain T> in the Hilbert space TL^. 

Proof. Let $ e P. The vector-valued distribution : (p Ll (xi) ■ ■ ■ ip ln {x n ) :\P is a re- 
striction of S l ^/8f l , I = Y2i=i^i> ^0 a partial diagonal with wave front set con- 
tained in {J xG M n i x } x V_ n . Since does not contain elements of the form 
[x\, k\, . . . , x n , k n ) with fej G V+, i = 1, . . . , n, the product is a well defined vector- 
valued distribution, and after smearing with some test function one obtains again 
a vector in V. ■ 



In particular, formula (22) makes sense everywhere, provided the expectation 
values of all time-ordered products of sub- Wick polynomials satisfy P4. More- 
over, every expansion into a sum of products of Wick polynomials by numerical 
distributions which satisfies ( |25| ) is of this form: 

P5. Causal Wick Expansion. 

T(L 1 (x 1 )---Z(x n ))= (^u,,T(^ l} (xx)---Z^(x n ))n u )x 

jl,— ,jn 

-.^(xt) ■ ■ '^"(in): 



jl! ■ ■ 'in! 



(26) 



4. Inductive Construction up to the Small Diagonal 

The properties defined in the previous section allow us to set up an inductive 
procedure in the spirit of Epstein and Glaser. We rely on a variation of their 
construction proposed by Stora [|4| . We start with a linear space W of Wick poly- 
nomials which contains all respective sub- Wick polynomials and want to define the 
time-ordered products T(Lx(xx) ■ ■ ■ £j(x n )) as a family of operator-valued distribu- 
tions which are multilinear in the entries £j G W and satisfy the properties Pl-5. 
We start the induction by setting T(l) = 1 and T(L) = L and assume that the 
time-ordered products for 1 < I < n — 1 factors have been constructed and satisfy 
all the defining properties. In a first step, we aim at constructing time-ordered 
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products of n factors on M n \A n where A n is the small diagonal submanifold of 
M n , i.e. the set of points (xi, . . . , x n ) with the property x± = X2 = • • • = x n . 

We use the space-time notion of causality in order to define a certain partition of 
unity for M n \A n : Let us denote by J the family of all non-empty proper subsets 
I of the index set {1, . . . , n} and define, accordingly, the sets Cj = {(xi, . . . , x n ) G 
M n I Xi J~(xj),i G I,j G I c } for any I € J. Note that the defining relation for 
the C/'s is related to causality on M and no£ on M n . It is fairly easy to show that 

Lemma 4.1. Let M be a globally hyperbolic space-time, then it holds 

|J C z = M n \A n . 
iej 

Proof. The inclusion U/Cj C A/ n \A„ is obvious. The opposite inclusion is proved 
as follows. Consider any set of points (xi , . . . , x n ) such that xi ^ Xj for some i ^ j, 
then the points Xi and Xj can be separated by a Cauchy surface £ as follows from 
the global hyperbolicity assumption. One may choose it as to contain none of the 
points Xk, k = 1, . ..,n. Hence, defining / = {k \ Xk G and noting that 

I G J we find {x\, . . . , x n ) G Ci. ■ 

We use the short hand notations 

T I (x I )=T([[L t (x l )), x 1 = {x ii iel). (27) 
iei 

The first step now is to set on any Cj 

Tj(x) =T\xi) T /c (x/c) , (28) 

as an operator-valued distribution since according to the induction hypothesis and 
the fact that / is proper this is a well defined operation on D(Ci). 

We now glue together all operators previously defined on different elements of 
the cover. For this we need to prove a sheaf consistency condition. Indeed, different 
C/'s overlap but due to the causality hypothesis P3 and the causal Wick expansion 
P5 valid for the lower order terms, the following property holds: 

Proposition 4.2. For any choice of h,l2 G J such that Ci 1 fl C/ 2 7^ we have 

Ti 1 Ic^nc^— Ti 2 \c It nc l2 j 

in the sense of operator-valued distributions over M n \A n . 

Proof. Let l\ , I2 G J and x — (x\ , . . . , x n ) G Ci 1 C)Ci 2 . Using the causality property 
P3 which is by assumption valid for time-ordered products of less than n factors 
we find, 

T^(x h ) = T^ixj^) T^{x IinI ,) , 

( 29) 

T'^xj;) = T^(z /fnj2 ) T^Hxj^) , 

and similarly for T l2 andT /2 . Now the terms T 11 0/2 andT /in/2 commute. Namely, 
they are based on mutually space- like points, thus using the Wick expansion for 
these terms this follows from local commutativity for the Wick polynomials of the 
free scalar field tp. Hence from definition (|^) we get on Ci 1 n C/ 2 

T Il = T Iinl2 T^ nl2 T Iinl2 T^ nI ^ , (Eqn.(||) + P3), 
= T /2 T /2 , (Eqn.@), 
= T l2 , (Eqn.@). 
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Let now {fi}iej be a locally finite smooth partition of unity of M n \A„ subor- 
dinate to {Ci}i£j. We formally define, following ( p7| ) and (pq), 

°T(£<i(a;i) . . . £„(:r„)) = ^ /j T/ . (30) 
Hence, we get our first crucial result, namely, 

Theorem 4.3. The symbols °T are well defined operator-valued distributions over 
M n \A n which satisfy the defining properties on T> C 7i w . 

Proof. We first prove that the definition does not depend on the choice of the 
partition of unity. Indeed, let {fj}iej be another such partition. Consider x G 
M n \A n , and let K. = {I e J \ x e C/}. Then there exists a neighbourhood V of 
x such that V C Pl/g/cC/, and supp(/j) and supp(/|) do not meet V for all I ^ K.. 
In this case 

^2ifi - f'i) T i\v= ^2(fi - f'i) T i\v ■ 

However, on V Ti is independent of the choice of I £ K.. Since J2ieic f 1 = 
S/eK f'i = ^ on we arrrve a ^ the conclusion. Furthermore, an inspection of 
the formula readily gives that the operator °T is defined on the domain X>, the mi- 
crofocal domain of definition of the Wick monomials, because of induction. Hence 
property PI. 

As far as the symmetry property P2 is concerned we just observe that the 
permuted distribution °T (x\, . . . ,x n ) = T{L^n\{x„(y\) . . . £^(71) (^(n))), has the 
expansion, 

°^ = £/7I7 = £/; (7) t; (/) , 

where we used the fact that the set J is invariant under permutations, but = 
Tj and {f^tn}ieJ is a partition of unity subordinate to {C/}/ e j, so symmetry 
follows from the result of the previous paragraph about the independence of T on 
the choice of the partition of unity. 

Causality P3 follows from an argument similar to the one used for the indepen- 
dence from the partition of unity. Indeed, taken any point x G M n \A n , as before 
x E V C DieicVi- From (§(]), 

T(x) = J2fi(x)Ti(x) . 
ieic 

Since Tj \y does not depend anymore on / G K, and 53zpx; fi = ^ over V hence 
°T(x) = Ti{x) which from ( |28| ) satisfies causality by definition. 

Now, we want to show that property P4 holds on M n \A n . It is sufficient to 
check that this property is satisfied for each Tj on Cj. We apply Wick Theorem to 
the components of the product in the definition (|2^) . It can be easily checked that 
the distributions ti in the Wick expansion of Tj are sums of terms of the following 
form 



flix^iXiy^Xio)- J[ U 2 (x i ,X j ) a ^ , 

(i,i)e/x/<= 



(31) 
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with dij S No, and where t 1 , t 1 are expectation values of lower order time-ordered 
products. 

The wave front set of ( |3l| ) is contained in the convex combination of the wave 
front sets of its factors. Hence it is given in terms of immersions of graphs with 
vertex sets I,I C , resp., and of a^j graphs with vertex sets {i,j}- All these im- 
mersions satisfy the condition in P4, the first two by assumption, and the last 
ones because of the definition of Cj and the properties of the wave front set of a 
Hadamard state. The union of these graphs is a graph with vertices {1, . . . , n}, and 
any convex combination of the components is given by an admissible immersion of 
this graph. 

Finally, property P5 follows from expression ( |28| ) by a straightforward applica- 
tion of the (generalized) Wick Theorem. ■ 



5. Steinmann Scaling Degree and the Extension of Distributions 

We now want to extend °T(Li(xi) . . .£j n (x n )) to the whole M n . As discussed 
before the problem can be reduced to the extension of the numerical time-ordered 
distributions °t(xi, . . . ,x n ) = (Cl u) ,°T(Li(xx) ■ ■ . & n (x n ))Q u ) . 

The extension can be performed in two steps. First °t is extended by continuity 
to the subspace of test-functions which vanish on A n up to a certain order, and 
then it is arbitrarily defined on a complementary subspace. It is this last step which 
corresponds to the method of counterterms in the classical procedure of perturbative 
renormalization. The extension of °t by continuity requires some topology on test- 
function space. The seminorms used by Epstein and Glaser in their paper are quite 
complicated, and their generalization to curved space-times appears to be rather 
involved. We found it preferable therefore to apply a different method already 
introduced by Steinmann |53f , namely the concept of scaling degree at a point of a 
distribution (see also fuH ). Its generalization to curved space-time is very similar 
to the concept of the scaling limit as introduced by Haag, Narnhofer and Stein [^oj 
and further developed by Fredenhagen and Haag p4l. A similar technique is used 




On Minkowski space, by translation invariance, the distribution is in terms of 
relative coordinates everywhere defined up to the origin, and there the concept of 
the scaling degree at a point leads to a rather smooth and economic method of 
renormalization, see e.g. flq] where the relation to differential renormalization is 
elaborated. On a curved space-time one needs the corresponding notion for a scaling 
degree with respect to the submanifold A n , and one also needs some uniformity 
of the singularity along A„ as well as control of the wave front sets during the 
extension process. 

Our strategy will be that, at first we introduce this improvement in the case of 
R™, then we discuss the case of manifolds. There we try to set up a procedure 
which allows to restrict the discussion to the pointwise case. 

5.1. The scaling degree. For simplicity, we work at first on M. d . Hence, consider 
a distribution t 6 D'(M d ). Let the action of the positive reals (dilations) be defined 
via the map 

A:R + x T)(R d ) — ► T>(R d ) 

(A,0) > (f> x = A~ d <f)(A~ 1 • ) , 
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and obtain, by pull-back, the map over distributions t £ D'(R d ) as, 

(A*t)^)=t x (<j>) = t(^) , 
where this operation in case t G Ll oc (R d ) is given by the explicit formula, 

t x ((j)) = ( t(Xx)<j)(x)d d x , V0 e D(R d ) . (32) 



The map A is clearly continuous w.r.t. the topology of T)(R d ) and we shall some- 
times use the previous formula ((H), by the usual abuse of notation, also in the 
general case. 

We say that t has scaling degree sd(t) = ui w.r.t. the origin in R d , if ui is the 
infimum of all uj' G R for which, 

UmA w 't A = 0, (33) 

holds in the sense of T>'(R d ). 

It should be clear from the definition that every distribution t G CD'(R d ) has a 
scaling degree to G [— oo,+oo[. If the distribution is not defined at the point we 
want to check, then the scaling degree might also be equal to +oo. We give some 
examples. 

Examples . 

1. Trivial example. Every (f> G £(R d ) has sd(^) < 0. 

2. Dirac measure. Let /i G £'(R d ) with n(cj>) = 0(0), G £(R d ), then sd(/i) = d. 

3. Feynman propagator. In the case of a free massive scalar field which is co- 
variant under translation on Minkowski space-time, the Feynman propagator 
can be written as 

/pip-x 
~ 2--^' 
p z — m z + le 

from which it is readily seen that sd(Ep) = d — 2. 

4. Homogeneous distributions. If t G D'(M. d ) is homogeneous of order a at the 
origin, i.e. t\ = \ a t, then sd(t) = —a. 

5. Infinite degree. The smooth function x — * exp(l/x), x G R+, is not defined 
at the origin and its scaling degree w.r.t. the origin is clearly infinite. 

As inferred from the 4th example, the scaling degree may be seen as a general- 
ization of the notion of the degree of homogeneity. Actually, our extension method 
is similar to the extension to all space of a homogeneous distribution as discussed in 
Hormander's book j}5| which on the other hand is also quite similar to the Epstein 
and Glaser procedure of distribution splitting [^0| . The fact that homogenous ex- 
tensions not always exist is the mathematical origin of the logarithmic corrections 
to scaling found in renormalization. Here, a discussion about space-time symme- 
tries and their implementation after renormalization is absent. It will be presented 
mQ. 

Lemma 5.1. The scaling degree obeys the following properties: 

(a) Let t G T)'(R d ) have sd(t) = lu at , then 

1. Let a G N d be any multiindex, then sd(d a t) < uj + \a\ . 

2. Let a G N d be any multiindex, then sd(x a t) < u> — \a\ . 

3. Let f G £(M rf ), then sd(ft) < sd(i) . 

(b) For U G D'(IR dl )> i = 1, 2 we have sd(h ®t 2 ) = sd(ti) + sd(i 2 ) . 
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Proof. The first two cases in (a) as well as (6) are straightforward. The third case 
in (a) follows from the fact that, by the Banach-Steinhaus principle, a convergent 
sequence of distributions is uniformly bounded. Hence, for every u' > sd(t) and 
every compact set K cK. d there is some polynomial P such that 

\X u 'tx(<f>)\ < sup \P(d)<f>(x)\ = H0IU.P . (34) 

Hence, for / 6 £(R d ), we have 

K/*)a(0)| = |*x(A0)|<A w '||/a0|Up ■ (35) 

The statement follows now from the boundedness of the sequence ||/a0||oo,p as 
A^O. ■ 

5.2. Extensions of distributions to a point. We now want to show how to 
extend a distribution t G D'(]R d \{0}) to all space by using the concept of the 
scaling degree. The scaling degree can easily be defined for such distributions by 
restricting the test functions appropriately. Equivalently, we can also, for each 
X 6 £(K d ) with ^ supp(x) look at the behaviour of the sequences x*Ai now 
considered as sequences in D'(R d ). 

There are three possible cases; when the scaling degree is +oo; in this case no 
extension to a distribution on K exists; when the scaling degree to is finite, but 
u> > d; then a finite dimensional set of extensions exists; or otherwise ui < d. We 
first study the third case. 

Theorem 5.2. Let to 6 CD'(R d \{0}) have scaling degree co < d w.r.t. the origin. 
There exists a unique t G D'(M. d ) with scaling degree lu such that t(<p) — to((f>), 

0e D(R d \{0}). 

Proof. The uniqueness is easy. Indeed, the difference among two possible extensions 
would be a distribution with support at {0}. By a well known structural Theorem 
of distribution theory this last is given by P(d)5 where P is a polynomial of degree 
deg(P) and 6 is Dirac measure at the origin. But this distribution has scaling 
degree equal to d + deg(P), hence a contradiction. 

Let us now consider a smooth function of compact support i9 such that $ = 1 in 
a neighbourhood of the origin. Set i3\(x) = ^(Ax), A £ K and 

f(«) = (l-tf 2 „)i , neN, 

where now is a sequence of distributions defined on the whole M. d . We wish 
to show that the sequence converges in the weak* topology of D'(M d ). Because of 
the sequential completeness of D'(R d ) it is sufficient to prove that it is a Cauchy 
sequence. Let <fi £ D(M. d ) and look at 

(t(»+ 1 )-tW)(0) = (^o)( 1 ? 2 „-^ 1 ) 

j (36) 

= 2-" d (# ) 2 -n(tf-tf 2 ) • 



According to Lemma 5J, (a)3., this sequence is majorized, for every u)' £ ]to,d[ by 
const. 2 n ' w ~ d \ hence it is summable as required. The limit 

t{4>) = lim t (n) {4>) , V</> G D(K d ) , 

n — >oo 

then defines an extension of to. It is obvious that the scaling degree of t is not 
smaller than u>. It remains to proof that it is not bigger than w. 
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Pick </> G and consider the following expression: 

t x {<t>) = lim A- d t o ((l-i? 2 n)0A-O ■ 

n — >oo 

Let R, e > be such that supp(0) C {x, \x\ < R} and -&{x) = 1 for \x\ < e. Then, 

(1 -d{2 n x))(t>{\- 1 x) = , 

whenever 2~ n e > Ai?. Let us choose n x G N such that 2~" A e > XR > 2- {nx+1 ^e. 
We have, 

oo 

(37) 

= £ (2"A)- d (t ) 2 -„((^- a )fc-„ A -i) . 

The set {($ — $2)^1 A* < const.} is bounded in D(R d \{0}). Hence for every 
to 1 > lu we find a constant c > such that 

|(to)a-»((0-02)02-»A-OI< C2 ""' 

Inserting this estimate back to Eqn. ( j37| ) we have 

00 <) — n\(d—u> ) 

I*a(«I < £ 2_n(rf " w,) = cA ' d 1 _ 2 -(^ ) 

n=n\ 

< C -r^\- d fH^V A<^' < c'A-"' 



1 - 2-( d -"' 

for some constant c' > 0. This proves the assertion. ■ 

We now deal with the extension procedure in case a distribution has a finite 
scaling degree lu > d. This extension procedure corresponds to renormalization in 
other schemes. To adhere more to the standard notation we introduce the degree of 
singularity p = lu — d. This is the analog of the degree of divergence of a Feynman 
diagram. 

Let D p (M. d ) be the set of all smooth functions of compact support which vanish 
of order p at the origin and let W be a projection from D(M. d ) onto D p (R d ). Since 
the orthogonal complement of Dp(R <1 ) consists of the derivatives of the (5-function 
up to order p, W is of the form 

W<j> = cj> - ™ a d a m , (38) 

\a\<p 

with Xo a being smooth functions of compact support such that c^ro^O) = 6g. 

Theorem 5.3. Let to G I>'(R d \{0}) have a finite scaling degree lu > d. Then there 
exist extensions t G D'(M. d ) of to with the same scaling degree, and, given W, they 
are uniquely determined by their values on the test functions fo Q . 

Proof. Any (j> G ©(R'') can be uniquely decomposed as <f> — <j>i + (f>2 where <\>\ = 
E|a|<p»o9 a 0(O) and (j) 2 G £> p (R d ). cf> 2 has the form, 

\a\ = [p] + l 
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with e 1>(R d ). We set 

(t,4>)= (x a to,i>a) + (t,<fix) ■ (39) 

l«l=[p]+l 



Since, by Lemma 5.1, x a t has scaling degree equal to p— [p] — l + d which is strictly 
smaller than d this term has a unique extension by Theorem |5.2| . 
We now prove that t has the same scaling degree as to- We write, 

tx(<f>) = (to o W) x {4>) - (d a 4>) (0)A- d -l Q l . 

\a\<p 

The second term clearly has scaling degree less or equal to p + d = u>. The first 
term can be written in the form 

(t )x(W(f>) + (t o )x((W0 x -i)x - W4>) . 

By assumption, the first term has scaling degree u>. To analyze the second term we 
write 

((W(j> x -,)x-Wct>){x)= ]T - A-l«lm Q (A.x)) (d a cf>) (0) . 

M<p 

(Note that (»«(•) - \-^tO a (X •)) 6 D p (M. d )). Using the identity, 
(tr Q (x) — \~' a 'tv a (\x)) = / —p~' a 'to a (px) dp, 



dp 

(pxk(dktv a )(px) - \a\m a (px)) ^i~' a ' _1 dp , 

(40) 

we get, after a moment of reflection for the exchange of the order between integra- 
tion and duality, that, 

(d a cf>) (0) f p- d -^-\h) x ^{{x k d k - \a\)tt> a )dp . 

The integrand can be estimated according to Lemma 5.1 . Indeed, for any u>' > lu 
we have, 



| (*o)x M -i ({xkd k - |a|)tt) Q ) | < const. (A 1 p)^ , 



and therefore 
-l 



which proves the assertion. 



A 



< const. A - "' — ■ — — - 



The expert reader can now proceed from this point to study the renormalizability 
of any theory which admits space-time translation covariance. The ambiguity of the 
extension is given by terms localized over the origin. The coefficients of these terms 
can be fixed by additional requirements, as customary in perturbative quantum 
field theory. We refer the reader to |ll[] for more details. 

During this process, one needs estimates on the scaling degrees of the arising 
distributions, corresponding to the power counting rules. In addition to Lemma 
5.1 estimates on scaling degrees of products of distributions (provided they exist) 
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are required. These can be obtained by explicit calculations (see e.g. the analo- 
gous estimates in Q). Much more elegant is a general method which exploits a 
microlocal version of the scaling degree. This technique is actually necessary if one 
wants to generalize the methods above to generic manifolds. We shall describe it 
in the next section. 

6. Surfaces of Uniform Singularity and the Microlocal Scaling 

Degree 

The generalization of the previous procedure to the case of submanifolds is what 
we really need in the treatment of perturbation theory on curved spaces. Indeed, 
the description given in Section 4 led to the notion of a scaling degree w.r.t. the 
small diagonal A n of the topological product M n . Here we classify the behaviour 
of distributions near some surface by a microlocal version of the scaling degree. 
We introduce two different notions. The first one, the (microlocal) scaling degree 
at a surface, involves only the surface under consideration, the second one, the 
transversal scaling degree, involves a fibration of the surface by transversal surfaces. 
The first notion behaves very nicely under tensor products and restrictions, whereas 
the second one admits an easy generalization of the extension procedure. As a 
matter of fact, the notions can be shown to be equivalent. 

6.1. Scaling degrees at submanifolds. Let M be a smooth paracompact mani- 
fold of dimension d and t be a distribution in D'(M). Let N C M be a submanifold 
such that the wave front set of t is orthogonal to the tangent bundle T3\f of N, i.e. 
for (x, k) S WF(i) with k 6 T*M, 

(fc,0=o, veer^. (4i) 

Under these circumstances, t can be restricted to a sufficiently small submanifold 
GcM which intersects 3\f in a single point x$ , such that the intersection of their 
tangent spaces at xq is trivial and their sum spans the whole tangent space (the 
submanifolds 6 and N are transversal, see e.g. [^7j, symbolically 6 rh N). This is 
due to the fact that WF(t) does not intersect the conormal bundle N*C = {(x,k) £ 
T*M\(k,£) = 0,V£ € T X G} of 6. Namely, for k e T* o M, (x ,k) £ WF(t) we have 
(k,£) = for £ e T Xo 7i, hence ^ for some £ € T Xo <5, thus (x ,k) (£ N*G. 

But WF(i) fl N*C is a closed conical set in TgM, hence its complement is an open 
conical neighbourhood of T* Q JA, in particular it contains a set Tfj M where Uq is 
an open neighbourhood of xq in C. By choosing C = Uq we arrive at the conclusion. 
So we proved, 

Lemma 6.1. Let t G D'(M) be a distribution on a smooth manifold M and let 3\T be 
a submanifold such that WF(t) _L T3sf. Then t can be restricted to every sufficiently 
small submanifold C such that IN" rh C . 

The singularity of t \e at xq may be classified by a covariant extension of the 
notion of the scaling degree, or better by a slight extension which uses microlocal 
analysis. For the economy of the presentation we first look at the concept of scaling 
degree at some surface N which reduces for each transversal surface C to the scaling 
degree at the intersection point. This last will just be a pointwise reduction of the 
general case we proceed to discuss right now. 
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Let U be a star-shaped neighbourhood of the zero section Z(T^M.) and consider 
a map a : U — > a(U) cNxM which is a diffeomorphism onto its range and such 
that the following properties hold true 
(i) a(x, 0) = (x,x), x G K; 
(m) a(TXn[/) cNxN; 
(m) a(sf,£) e {x} x M, z e K, £ 6 T X M; 
(w) dea(af,-)^ = o=idT a> M- 

A concrete example of such a map a can be defined, whenever we consider the 
manifold M endowed with a (semi-)Riemannian metric, in terms of the exponential 
map, namely, a(x, £) = (x, exp x £), provided the submanifold N is totally geodesic, 
as will be the case in our applications. In the general case, we shall call the set of 
all such maps by Z. 

Let a G Z and set t a = (l<g>t)oa on T>'(U) and t%(x,£) = t a (x, A£), < A < 1. 
Here, (1 <g> t, <t> ® tp) = J <j) ■ (t,ip) for test-densities <f> G and ^ 6 T>i(M). 

Since U is starshaped, A _1 C/ DC/ for < A < 1, hence can be considered as a 
distribution on T>i(U). 

As a preliminary step we have the following 



Proposition 6.2. For any t £ D'(M) which satisfies the hypothesis of Lemma 6.1, 
there exists a closed conic set V <zT*U such that 

{i) r _L T(TNn U); 
(m) WF{tf) C r. 

Proof. Since a maps TNCiU into NxN, its derivative a* : TC7 — > T(N x M) maps 
T(T^n U) into T(N x N). But WF(t) 1 TN implies WF(l®t) _L T(N x N), hence 

WF(f) = a*WF((l (8 t)r a(t0 ) -L a^T^N x N)) = T(TNn £7) . 

Now, 

WF(^) = {(*:,£; fc) e T^^K^A^fe) £ WF(f)} . 

Here, we identified the cotangent spaces at the points (x, £) and (a;,A£) by the 
isomorphism induced by the diffeomorphism U — > AC/, (x, £') — > (x, A£'), £' G T^M. 
Now, let £ G T X N and 77 G We may identify n with a vector in Tr x xgfli 

and observe that it is orthogonal to WF(t") and hence also to WF(t°). We now 
set r = Uo<A<iWF(i°), where the closure is performed within T*U . It remains to 
prove (1). 

Let (x,£„;fc„) G WF(t^) be a convergent sequence in T*U with limit fx, 
k ^ and £ G T^N. There is a corresponding sequence (x,A„£ n ;fc„) G WF(t a ). 
Let A G [0,1] be a limit point of the bounded sequence {A„} ng N- Then, there is 
a subsequence converging to (x, A£;fc) G WF(<"). But A£ G T X N, hence we have 
k _L T( x ^)(TN). If we again identify the tangent spaces at (a;, A£) and (x,£) we 
obtain the desired result. ■ 

Choosing first any map a G Z we are ready for the following: 



Definition 6.3. A distribution t G D'(M) has the microlocal scaling degree u> at a 
sm 

if, 



submanifold N w.r.t. a closed conical set Tq C N*N, symbolically ui = /UsdoJ 1 (t, a) 



(i) i/jere exists a closed conic set V C T* (T^M) tuit/i i/ie properties stated in 



Proposition 6.i, with the first one replaced by T\tnC a*(Z(T*J{) X Fq), 
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(m) uj is the infimum of all those lu' for which, 



limA a; i? = 0, (42) 

in the sense of the Hormander topology on T>' T (Tj^M) . 

Now, depending on the position of Tq one can give different refined versions of 
the microlocal scaling degree. Indeed, when the inclusion in N*N is proper we 
speak of the strict microlocal scaling degree. When To = N*N, we call it simply 
the scaling degree at the submanifold, symbolically sd^(i). Moreover, when To = 
we speak of the smooth scaling degree. 

The definition seems to depend on the choice of the map a £ Z. In our concrete 
case we could make use of the metric to choose a canonical diffcomorphism a 
in terms of the exponential map, but for reasons which will become clear in the 
following it is helpful to proof its independence. 

Before coming to that point we show an example for the computation of the 
scaling degree which is relevant for the physical discussion, namely, the generaliza- 



tion of the example (3) in subSection 5.1, the Feynman propagator Ep, to curved 



space-time. The Feynman propagator is considered as a distribution on M x M , 
and we are interested in the microlocal scaling degree at the diagonal A2 C M x M. 
Indeed, the wave front set of Ep is orthogonal to the tangent bundle of the diagonal. 
We choose a : Ta 2 M 2 — ► A 2 x M 2 ~ M x M 2 as a(x, £1, £2) = (x, exp^. £1, exp x £2) 
and obtain as on Minkowski space sdA 2 (EF) = d — 2. A similar result holds for the 
2-point Wightman distribution u>2- So, 

Lemma 6.4. The microlocal scaling degrees of the 2-point Wightman distribution 
lo 2 at the diagonal A 2 w.r.t. T a = {(x, k; x, — k)\x £ M, k £ dV+, k ^ 0} is given by 
/isd^° 2 (w2) — d — 2, where d is the dimension of the space-time. 

6.2. Invariance and properties for the scaling degrees. Let us then choose 
two maps ai, (X2 € Z and state the following: 

Proposition 6.5. Let t £ D'(M). Let u>i = fisd^(t, Oi), i — 1,2 be the microlocal 
scaling degrees w.r.t. N and Fq resp. for the two arbitrarily chosen maps. Then 

Wi = Cl>2 ■ 



Proof. It is simple to check that 
where (3\{x,£) = X~ x P(x, \£) and (5 — a^ 1 o c*2- Now, assume u\ is the scaling 



CW=C(^^) . V0£ D^M) , (43) 



degree for t w.r.t. N and To associated with a\. We should proof that Eqn. ( |42| ) 
for a2 converges in the sense of D'(T-xl\t) as well, at the same rate as A J. 0. 

The convergence in the sense of distributions is simple. Indeed, if supp(<^>) C K, 
K a compact subset, then there exists a Ao such that /3a (supp (</>)) C K for all 
A < A . Hence it suffices by Banach-Steinhaus principle to proof that the family 
{<f> o Z?^ 1 ^ < A < Ao} is bounded, uniformly in A, in Di(K) w.r.t. the family 
of continuous seminorms which gives the appropriate Frechet topology. This check 
proceeds easily from the chain rule and the verifiable fact that the only contribution 
comes from the 0-th and 1-st order derivatives w.r.t. Z?^ 1 . In the limit they are 
the only terms which survive giving resp. the identity map on the bundle and the 
derivative of the identity map. Hence, we get the same rate of convergence as far 
as plain distribution convergence is concerned. 
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A little bit trickier is the convergence in the sense of seminorms for M2 (b). 
Starting again from Eqn. via the multiplication of a smooth test function of 

compact support ip such that ip = 1 on a small neighbourhood of supp (</>), we have 
that ipt is of compact support and then, by a partition of unity with functions with 
support on charts, that we are working on M. 5 x M. d , where S is the dimension of the 
submanifold N. Now, let us multiply the test function 4> with the term exp(i(fc, • }), 
and use inverse Fourier transform to get, 

= / ^( P )h(p,k;(3 x )d d+5 p 7 

where, 

l4(p,ki0x) = fe«Me>*>-<t*»m d d+5 ^ 

where in all these expressions the coordinates £ are the local coordinates of Tj^M 
and k and p are their dual coordinates. 

We use the idea of the proof for the stationary phase Theorem, see for instance 
Theorem 7.7.1 in Hormander books fl35| . Because of f3\ — > id for A — > the 
oscillatory integral 1$ falls off rapidly outside of any conical neighbourhood of the 
diagonal p = k in R d+A x R d+d , uniformly for A sufficiently small, i.e. for every 
e > there exists a Ao > such that for every N £ N 

sup sup (1 + \p\ + \k\) N \l4,(p, k, p x )\ < oo . (44) 

0<A<A |p— /s|>e|fe| 

Now let T c M. d+S be a closed cone such that 

sup(l + \p\) N \W?(p)\\ u - , A - , (45) 
c 

for all closed cones C with C l~l V = and all N £ N. 

We now want to show that the same property holds for 4>t" 2 . So let C be a 
closed cone such that the closed cone 

C = {p£ R d+S , \p-k\< e\k\ for some k £ C} (46) 

does not intersect T. Then we split the region of integration over p into the parts 
|p — k\ < e\k\ and the rest. In the first region we can estimate k by p and use the 
fast decay of (1 + I'l/jt'? 1 (p)\ within C and the polynomial boundedness of 7^; 
in the second region the polynomial boundedness of (1 + IpD^I?/^" 1 (p)\ and the fast 
decay of 1$ (p|). This proves the desired estimate for (f>t^ 2 . ■ 

The microfocal scaling degree /^sd has similar properties as the scaling degree, 
as described in Lemma |5.l| . In addition, one finds the following two properties: 

Lemma 6.6. Let ti,ta £ D'(3VC) with nsd u>\ resp. u>2 at N C M, w.r.t. Tg resp. 
Tg and such that Z(N*yi) £ (Tj + Iq). Then the pointwise product t±t2 exists in a 
small neighbourhood of N and has the microlocal scaling degree lu < u>i + u>2 at N 

tu.r.t r = rSurgu(rJ + r§). 

Proof. By assumption, the wave front sets of i" A and t% ^-i some a £ Z, on a suf- 
ficiently small neighbourhood of N satisfy the condition (WF(<" A ) + WF(i" j)) R 
Z(T*(TyjM)) = 0, hence their product exists there by Ml. Because of the se- 
quential continuity of the products in the Hormander topology M2, the microlocal 
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scaling degree is given by the sum w.r.t. the stated conic region as follows from 
Ml and does not depend on the choice of the map a. ■ 

The following nice property follows from the sequential continuity of the restric- 
tion operator to submanifolds M3: 

Lemma 6.7. Let Ni be a submanifold o/N, and let t G D'(M) have the microfocal 
scaling degree ujatJ4 w.r.t. Tq. Then the microfocal scaling degree oft at Ni w.r.t. 
the restriction Ti o/Tq to Ni is less or equal to to, /isd^^) < ^sd^(i). 

A last word is devoted to a pointlike trivialization of the above procedure. This 
case can be derived straightforwardly by considering N = {p}, where pglisa 
generic point and thought of as (a rather singular case of) a submanifold. 

The translation to this simpler case is done via the following correspondence 
between geometrical and analytical quantities: 

U G TjstM — ► Up G T P M , 
a : U — » N x M — > a : U p — > M , (i) and (iv) valid, 
t a = (1 ® t) o a — > t a = t o a , 

r c t*u , r i T(TNnu) — > r c t*u p , r i tu p , 

T Q G N*N > T p G r p *M . 

6.3. Transversal scaling degree. Instead of blowing up distributions on M to 
distributions on N x M in the definition of the scaling degree, one could also use 
a fixed fibration of a neighbourhood of N in M by transversal surfaces. For this 
purpose we decompose T^M into complementary subbundles, Tj^M = TN+C. The 
map ac '■= ^2°o\cnV i with the projection 7T2 : 3\TxM — > M onto the second factor 
and with V being a sufficiently small neighbourhood of the zero section Z(TyjJA), is 
then a diffcomorphism onto some neighbourhood of N. The images of the fibers of 
C are transversal surfaces. The transversally (w.r.t. a and C) scaled distribution 
is then defined by 

t\,_L(x, n) =to a c {x, Xn) , (x, n) € C . (47) 

The transversal microfocal scaling degree may then be defined as the infimum of 
all uj G K such that the sequence X"t\,± converges to zero within D' Tc (C), with a 
closed conical set Tc G T*C with Tc\z(c) = a c(^o)- Fortunately, it turns out that 
this new concept of a scaling degree at a surface coincides with the old one. Thus, 
in particular, the transversal scaling degree does not depend on the choice of the 
fibration. 

Proposition 6.8. The transversal microfocal scaling degree defined above coincides 
with the microfocal scaling degree defined in (pl[). 

Proof. We may restrict ourselves to a sufficiently small neighbourhood of a point 
at the surface INT. In suitable coordinates, M is a subset of R s x R d_<s where the first 
factor corresponds to N and the second factor to the transversal surfaces. T^M is 
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a subset of 1* x I s x R d ~ s with the first factor corresponding to N, the second to 
the tangent spaces of N and the third one to the fibers of the transversal bundle 
C. The map a may be chosen as 

a(x,C,v) = (x,x + £,rj) . (48) 

Then etc becomes the identity. The distribution t may be replaced by a distribution 
with compact support, and the factor 1 in the blow up of t may be replaced by a 
test function \ G D(N). For the Fourier transforms we then obtain 

f x (p,q,k) = \- d x(p-\- 1 q)t(\- 1 q,\- 1 k) , (49) 

and 

t^l(p,k) = X d - 5 t(p,\- l k) . (50) 

Furthermore, using a corresponding trivialization of the respective cotangent bun- 
dles, we may identify To — Tc with {0} x K, where K is a closed cone in R d_5 , 
considered as the transversal part of the cotangent space, and V with {0} x {0} x K. 
The convergence of t\ may be discussed in terms of seminorms of the form J v (l + 
\p\) \t\\ with, respectively, conical neighbourhoods V of T and some N G (— N) 
and closed conical sets V in the complement of T and all N G N. Since \ is strongly 
decreasing, these seminorms of t\ can be estimated in terms of the corresponding 
seminorms of t^±. I 

6.4. Extension of distributions to surfaces. We now want to apply these con- 
cepts to the extension problem of distributions t G D'(M\N). The wavefront set 
of the extension shall be orthogonal to TJvf, hence a necessary condition is that 
this holds true for the closure of WF(i) within T*M. We extend the notion of 
the microlocal scaling degree to such distributions in an analogous way as in the 
extension problem to a single point. Namely, for an arbitrary function \ £ £(M) 
with suppxHN = 0, (1 ® x) ° ct ■ can be considered as a distribution on U C T^M. 
The microlocal scaling degree of t at Jvf is then defined in terms of all sequences so 
obtained. 

We choose a fibration of a neighbourhood of N by transversal surfaces. It is easy 
to see that if t has a scaling degree to at N, than its restriction to a transversal 
surface has a scaling degree at the point of intersection with 3\f which is less or 
equal to tu. We therefore obtain the corresponding extension theorem. 

Theorem 6.9. Let N be a submanifold of the manifold M, and let to G D'(M\?sf). 
(i) Ifsdj~i(to) < codim(N) there exists a unique distribution t G D'(M) extending 

to with sdyi(t) = sd>r(to) 
(ii) If codim(3\f) < sdH(^o) < 00 there exist extensions t G D'(M) with sd^(t) = 
sdjvf(to). They are uniquely characterized by their values on some closed sub- 
space o/D(M) which is complementary to the space of all test functions which 
vanish on N up to order sd>f(to) — codim(!N}. 

Proof. According to Theorems |5.2| and |5.3| there exist extensions of the restrictions 
of to to every transversal surface with the same scaling degree. We have to show that 
they are restrictions of a unique distribution t on M with the same microlocal scaling 
degree. We first fix some normal fibration as described above and consider to as a 
distribution on ?7\{0} with a neighbourhood U of the zero section of C. We perform 
the construction of t at all fibers, by choosing a smooth function § G £(C) which 
is equal to 1 in a neighbourhood of the zero section and whose restrictions to every 
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fiber have compact support. Moreover, we choose smooth functions to p £ £(U), also 
with compact support on each fiber, which satisfy the condition djtt)p(x,0) = 5^ 
where £ denotes the variable in the fiber over x £ N. We may take top — tv^/(3\ 
with some function tu which is identical to 1 in a neighbourhood of the zero section. 

We set p = sdj~i(t ) - codim(3\f), # A -i(a;,£) = ?9(a;,A _1 ^) and W(j){x,^) = 
0( :z: >O-E|/3|< (9 tt '/3(2 ; ^)'9^(a;,O). Let T c = WF(t )\JN*N. We already know that 
the sequence t n = t a {l — ^ 2 ™) W converges on every fiber, and it is easy to see that 
it converges weakly in D'(U). We now want to show that the wave front set of t is 
perpendicular to TN. For this purpose we show that the above sequence converges 
even in D' Tc (U), i.e. that for every pseudodifferential operator A whose wave front 
set does not intersect Tc (hence At n is smooth) the sequence At n converges in the 
sense of smooth functions. For a pseudodifferential operator with smooth kernel 
the argument is essentially the same as for the weak convergence, hence we may 
restrict ourselves to pseudodifferential operators whose kernels have support in a 
sufficiently small neighbourhood of the diagonal of U x U (only here singularities 
may occur). 

According to the discussion of condition M2(b), we may equivalently look at the 
Fourier transform of x^o(^2 m — ^2 m +0 where x is a test function with sufficiently 
small support which does not vanish at some point x £ 3\f. Introducing suitable 
local coordinates in a neighbourhood of x , we find 

[tfoOV. - tf 2m+ 0] A (P, k) = 2- m ^ [(xto) 2 -m±(l? - 2 )] A (p, 2- m k) , 

where the symbol [ • ] A means Fourier transform, the p's are the dual coordinates 
w.r.t. points x £ INT and similarly the k's are dual coordinates w.r.t. the points 
£ £ C x of the fibers C x of C and finally 5 is the dimension of N. 

By the assumption on the scaling degree of to at Jsf we know that for all e > 0, 
JVgN and ui > sdx(to) there exists some c > such that 

|[(xfo) 2 -m,±(l? - #2T(P, k)\ < c 2 m -(l + |p| + Ifcl)^ , (51) 

for all (p,k) with \p\ > e\k\ (i.e. outside of a certain conical neighbourhood of the 
normal bundle {(p, k),p = 0}). Therefore, the sequence 

sup |[xto(0 2 .» -i? 2m+ i)] A (p,fe)|(l + H + |fe|) JV , (52) 

\p\>e\k\ 

is summable for oj < (d — S). If sdj^(io) < (d — S) such an u> exists, and we conclude, 
that the sequence t (l — ^2 m ) converges in D' Tc (M). 

In case of sdn^o) > (d — 5) we have to apply the ^-operation defined above, 
which, for \ with sufficiently small support, reduces to a subtraction of the Taylor 
series up to order [p] . We obtain from the integral formula for the remainder in the 
Taylor expansion 

{ x t (d 2m -tf 2m+1 )oW} A (p,k) = 

f E ^o(^-^ + oer(P,^ (1 ~ M)[P ' (M + 1) ^. 

Jo l/3|=W+i P] 

Since X^o^ 3 has scaling degree sd^-(to) — [p] — 1 < (d — 5), we may use the same 
estimate as before and find that the sequence to(l — "&2 m ) ° W converges in D' r (M). 

It remains to prove the stability of the scaling degree under the extension pro- 
cedure. The argument is a straightforward combination of the techniques in the 
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corresponding proofs in Section 5 and the arguments above and is therefore omitted. 



7. Extension to the Diagonal and Renormalization 

We come to the main point, namely to prove that the inductive analysis of 
Section 4 closes when supplemented by the information about the scaling degrees 
and gives well defined operator- valued distributions T n all over the space M n . This 
process is what it is usually called "renormalization." 

Toward this goal we must show that the scaling degree for the distributions of 
the nth order of the induction can be estimated in terms of those of lower orders. 
As explained before it is sufficient to do it for the numerical distributions 

°t(xi, ...,x n )= w(°T(xx, . . .,#„)) ! (53) 

where to is our reference Hadamard state. According to ( fjl] ) °t is a finite sum of 
terms 

fi{x)^{xi)^\xic) n wafo.ay) *' , (54) 

with nonnegative integers ay and a smooth partition of unity of M n \A n fj G 
£(M n \A n ) with supp// C Cj. M n inherits a natural metric from M, and all partial 
diagonals Aj are totally geodesic submanifolds. The map a : TM n — > M n x M n 
may therefore be defined by a(x,£) = (x, exp x £). Then all restrictions of a to 
partial subdiagonals, a>i = o.\t a M n satify the conditions before Proposition 6.4. 
We can choose the functions fj with smooth scaling degree at the small diagonal. 
We then may consider all factors in (|54|) as distributions on M n . According to 



Lemma 6.7, their microlocal scaling degrees with respect to A„ are bounded from 
above by their microlocal scaling degrees with respect to the respective partial 
diagonals. Moreover the convex combinations of the respective conical subsets of 
the con orm al bundle of the small diagonal do not meet the zero section. Hence, by 



Lemma 6.6, the scaling degree of the distribution in ( |54| ) is bounded by 

sd A , (ti) + sd Af c (tic ) + ^2 a v ( d ~ 2 ) • ( 55 ) 



us 



From this we get the formula 



n n , _ _ 

sd A „ ( W (T(: J] <P li (Xi) :))) < ^ k— . (56) 



i=i 



We thus obtain 



Theorem 7.1 (Main Theorem). All polynomially interacting quantum field theo- 
ries based on the the scalar field on d > 2 dimensional globally hyperbolic space- 
times follow the same short- distance perturbative classification as on their respective 
Minkowskian cases. 

We recall that for simplicity we have considered only pure monomials as inter- 
acting terms, i.e. without derivatives, multiple interacting fields and so forth. The 
general case however can be derived straightforwardly from our construction and 
we leave the task to the reader. 

We close this section by specifying the choices which have to be made in the 
abstract geometrical setting of Section 6. First we choose the normal fibration of 
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M n in a neighbourhood of A n which was used in the proof of the Theorem 6.9 for 
the construction of the extension. Let iVA„ be the orthogonal complement of TA n 
in T A „M n w.r.t. the metric in M° 

NA n = {(*,&,... ,£n) eTA„M",^ei = 0} . (57) 

Then 7:2 ° c4jva„ describes the desired fibration, 

7T 2 oa(i,(i, ... = (exp^^i, . . . ,exp x £„) . (58) 

x may be considered as the center of mass of the points Xi = exp x ^ , and the tangent 
vectors £j with the constraint = pl & y the role of relative coordinates. We 
further choose a smooth function tn on TM which is equal to 1 on a neighbourhood 
of the zero section and has compact support on each fiber and use the function 

tt) n (af,Ci,... ,£„) = n tt, ( a; '&) > ( 59 ) 

in the extension to the diagonal A„. 

By these conventions we get a reference definition of time-ordered products for 
all Wick products which involves only the chosen Hadamard state and the function 
ro. The algebra of interacting fields can then be defined by choosing the coefficients 
in the Lagrangian. 

8. On the Definition of the Net of Local Algebras of Observables 

In the preceding Section we finished the construction of time-ordered products 
of Wick polynomials of the free fields. We now want to show that this already 
gives the full net of local algebras of observables (within perturbation theory). An 
"adiabatic limit," whatever this might mean on a curved space-time, is not required. 
Actually, these observations are not completely new. In |5^] it was already observed 
that the local S'-matrices of the Stiickelberg-Bogoliubov-Epstein-Glaser approach 
give a local net of observables. 

Let W be the set of Wick polynomials with coefficients from D(M). So every 
A G W is an operator valued distribution with compact support which is relatively 
local to the free field. Our starting relation is the causal factorization (^) 

S(A + B + C) = S{A + B)S(B)- 1 S{B + C) , (60) 

for A, B, C G W, whenever supp(^4) is later than supp(C). 

Now let L be our interaction Lagrangian. Then gL G W for g G D(M). We may 
define observables with respect to the interaction gH by Bogoliubov formula 

S gL (A)^S(gL)- 1 S(gL + A). (61) 

We now show that the interacting observables depend only locally on the interac- 
tion. More precisely, we have the following 

Proposition 8.1. Let be a causally closed region, and let the test functions g 
and g' coincide on some neighbourhood of 0. Then there exists a unitary V such 
that for all A e W with supp(A) C 

VSgc^V- 1 =S g ,c(A) . (62) 

Proof. We may split g' — g = a + b where a does not intersect the past of and b 
not the future. Then supp(a£) is later than supp(j4), hence from ([30|) we find 

S(g'L + A) = S(g>Z)S((g + 6)iL)- 1 5(( 5 + b)L + A) , (63) 
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thus Sg'a(A) — Sfg+fyziA). Moreover, supp(A) is later than supp(6£), hence 

S((g + b)L + A) = S(gL + A)S(gL)- 1 S((g + b)L) . (64) 

Hence we obtain ( |62| ) with V = Sg^^L)^ 1 . ■ 

We conclude that the algebra 2l s £,(0) which is generated by S g n(A), supp(A) C 
with g = 1 in a neighbourhood of 0, is up to unitary equivalence uniquely fixed 
byL. 

We may formalize the construction in the following way. Let 0(0) for a causally 
closed compact region be the set of test functions which equal unity on a neigh- 
bourhood of 0. Consider for g,g' £ 6(0) the set V sff /(0) of unitaries V which 
satisfy the intertwining relation 

VS gL {A) = S a , L {A)V , supp(A) c . (65) 

The algebra of observables 21^(0) can now be defined as the algebra of covariantly 
constant sections of the bundle 

(J i9}x%c(0). (66) 

see(o) 

Here, a section A = (^4 g ) ff ee(0) is called covariantly constant if 

VA g =A g ,V, WeV gg ,(0). (67) 

21^(0) contains for example the elements Sc(A) = (S g &(A)) ge ®r y To complete 
the construction of the net of local algebras of observables we have to fix the imbed- 
dings «o 2 Oi : 2l,c(0i) — > 2l,c(02) for Oi C ©2- But this structure is inherited from 
the fibers and may be defined by the restriction of the section from 0(Oi) to 0(02)- 

9. Summary and Outlook 

We have proven that renormalization on curved backgrounds can be done in 
close analogy to renormalization on Minkowski space, and that the removal of sin- 
gularities follows the well known power counting rules. This result was expected 
since the ultraviolet behaviour on smooth manifolds should be essentially identical 
to that on Minkowski space. But we had to overcome two major obstacles: On a 
generic space-time there is no reason to expect a decent infrared behaviour, hence 
we had to use a method which decouples completely the short distance from the 
long distance problem; the other source of problem was the absence of transla- 
tion invariance which, on the technical side makes obsolete the usual momentum 
space methods, and on the side of physics, forbids to base the construction on a 
distinguished vacuum state and on the notion of particles. 

We solved these problems by basing the construction on the local S'-matrices 
of Stiickelberg and Bogoliubov, by invoking the general ideas of algebraic quan- 
tum field theory and by replacing translation invariance by smoothness properties, 
making extensive use of techniques and concepts from microlocal analysis. 

Besides the solution of the problem to which this paper is addressed, we solved 
several other problems which might be of independent interest. First we have 
found a new construction of Wick polynomials on a domain which depends only 
on the representation associated to a fixed Hadamard state but not on the state 
itself. Since according to Verch p7| , all representations induced by Hadamard 
states are locally quasiequivalent this amounts to an algebraic construction of Wick 
polynomials. 
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Second, we gave a perturbative construction of algebraic quantum field theory. 
In particular, we proved that the theory (in the algebraic sense), is completely fixed 
if it is known locally. Actually this holds independently of perturbation theory and 
might be a hint for a construction (in the sense of constructive quantum field theory) 
of asymptotically free theories. There the construction in small volumes seems to be 
possible |Q but the infrared problem poses, at present, unsormountable difficulties. 
The message of this paper is that the construction of the algebra of observables is 
nevertheless possible once the algebra of observables for small space-time regions 
have been constructed. The long distance behaviour of such a theory would still 
require an extra investigation, but it would be the behaviour of an existing theory, 
quite similar to the computation of spectra of Hamiltonians which have been shown 
to be self-adjoint operators. 

On the technical side we had to study the extension problem for distributions 
which are defined on the complement of some submanifold. This seems to be 
a natural mathematical question, and in a simple case it was treated in [^3| by 
similar methods. What seems to be new is our concept of the (microlocal) scaling 
degree at a surface which combines the condition of smoothness along a surface 
with a classification of the singularity in a transversal direction. 

The main open point in this paper is the fixing of the finite renormalizations. 
One expects that they can be chosen in terms of local functions of the metric, but 
a precise formulation meets a lot of problems. A similar problem was studied (and 
partially solved) in the definition of the expectation value of a renormalized energy 
momentum tensor of free fields by R. Wald |6l| . We hope to return to this problem 
in a future publication pp| . 
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